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Abstract 

We study meson-meson interactions using an extended q^q^i^g) basis that allows calculating cou- 
pling of an ordinary meson-meson system to a hybrid-hybrid one. We use a potential model matrix 
in this extended basis which at quark level is known to provide a good fit to numerical simulations of 
a q^(f system in pure gluonic theory for static quarks in a selection of geometries. We use a combi- 
nation of resonating group method formalism and Born approximation to include the quark motion 
using wave functions of a qq potential within a cluster. This potential is taken to be quadratic for 
ground states and has an additional smeared - (Gaussian) for the matrix elements between hybrid 
mesons. For the parameters of this potential, we use values chosen to 1) minimize the error resulting 
from our use of a quadratic potential and 2) best fit the lattice data for differences of Eg and n„ 
configurations of the gluonic field between a quark and an antiquark. At the quark (static) level, 
including the gluonic excitations was noted to partially replace the need for introducing many-body 
terms in a multi-quark potential. We study how successful such a replacement is at the (dynamical) 
hadronic level of relevance to actual hard experiments. Thus we study effects of both gluonic exci- 
tations and many-body terms on mesonic transition amplitudes and the energy shifts resulting from 
the second order perturbation theory (i.e. from the respective hadron loops). The study suggests 
introducing both energy and orbital excitations in wave functions of scalar mesons that are modelled 
as meson-meson molecules or are supposed to have a meson-meson component in their wave functions. 



1. INTRODUCTION 

Given the availability of both the numbers generated by lattice simulations of quantum chromodynam- 
ics and continuum models of the hadronic systems, an effective use of the numbers could be to improve the 
models through constraints of getting a least chisquare difference with the numbers for the correspond- 
ing discrete quarks and antiquarks configurations. Such lattice-improved models can then be reasonably 
used for all spatial configurations to eventually give dynamical predictions for experimentally measurable 
quantities like meson masses, meson-meson bindings and cross-sections and shifts (polarization potentials) 
to meson masses arising through meson-meson loops, etc. For one pair of quark and antiquark, a well 
established such use of lattice results is substituting in a Schrodinger equation a Coulombic-plus-linear 
quark-antiquark potential supported by lattice QCD calculations (see ref. [H El E] and others) for the 
ground state of the gluonic field between a quark and antiquark. Now that lattice results for excited 
state of the gluonic field are also available for years, even some dynamical uses of excited state gluonic 
field potentials have been worked out [H O [B] . 

Uses worked out by others are either limited to numerical calculations without an explicitly written 
excited state gluonic field potential or the potential used originates from flux tube [7] or string models [3 

Each of these approaches has its usefulness. What we add to this series of works is ourselves writing 
an analytical quadratic plus exponentially falling expression for the excited state gluonic field potential 
between a quark and antiquark and fit its parameters to the lattice data for the excited state gluonic 
field values available for discrete quark antiquark separations in [lOj . This is reported in our previous 
work as well. This work of us actually suggests and evaluates few other expressions for the excited 
state gluonic field potential as well. But the dynamical applications in it are present for a system whose 
valence quark contents are limited to one quark and one antiquark. 
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Through the present paper, we extend work on the dynamical imphcations of gluonic excitations to 
multiquarks that can be composed to more than one hadronic subclustcrs. For this extension, we need 
quark-level potentials that can model the more complicated gluonic field of this multiquark system. For 
this we combine the modelling of the spatial distribution of this gluonic field reported in ref. 11 (along 
with its fits to the continuum limits of the corresponding lattice simulations) with a realistic three colour 
structure. We have to do this combination because using all three colours the direct lattice simulations 
of all the Wilson loops relevant to a two quark two antiquark system are perhaps limited to ref. |12j and 
those mentioned in refs.[9-ll] of ref. [T3]- These works use a basis for a q^if system that is truncated 
to the ground state of the giounic field. In comparison ref. |11| extends the basis to include the gluonic 
excitations and its description of the spatial distribution of the q^cf giounic field is more complete. But 
the lattice simulations in ref. ,11] were carried out in a two-colour approximation to save computer time. 

Our purpose here is to take advantage of the relatively complete basis and spatial distribution models 
of ref. TT] but using all three colours in the quark-level potential we use for our calculations of the 
dynamical hadron-level implications for a meson-meson system. In changing a two colour based model to 
a full three colour one, the number of colour basis states that interact (for any the spatial configuration) 
remain the same (see eqs. A.l and A. 3 of ref. TJ], along with eqs. 5.1 to 5.4 and fig. 5.1 of ref. [TS]) 
and we had to essentially only replace some colour overlap factors with proper S'[/(3)c values as elements 
of matrices of the same order. But the inter-quark elementary potential of ref. [llj turned out to be 
problematic for including the quark motion for the hadron level implications and thus we had to replace 
their numerically fitted 0.562 + 0.0696r — "'^^^ — '^■'^2^ ground state quark antiquark potential by a 
constant plus quadratic confining potential term and the additional - — -f- 5^983 potential for the 
gluonic excitation by one of the form Aex-p{—Br^). As written in ref. [6 , only ^exp(— i?r^) can be 
used in solving the integrals of our present work analytically although few other forms for the gluonic 
potential are also suggested in ref. [6 with less a-s compared to Aexp—Br^. Only then we are able 
to perform a full meson level dynamical calculations for transition amplitudes from one set of quark- 
antiquark clusterings (mesons) to the other. Using these amplitudes we are also able to study certain 
properties of the polarization potentials for a meson-meson system. 

We had to use simpler interquark potentials that can be symbolically integrated at a later stage, after 
necessary multiplications by wave functions of qqqq positions, to complete the adiabatic-approximation- 
based treatment of a qqqq system mentioned below. Being not limited by such demands of later inte- 
grations, the form and then parameters values of the continuum qqqq model proposed in ref. [llj were 
chosen to simply minimize 

N{G) 

xi = ]^ E(^^-*^^)VAi?., (1) 

where N(G) is the number of data points for geometry G. The qqqq geometries in ref. [IT numerically 
worked on were (quarks at the corners of) squares, rectangles, tetrahedra and some other less symmetric 
geometries Q,N,P and L (linear). For each data point i , the lattice energy Ei was extracted by solving 
the following eigenvalue equation (see eq. 2 of ref. [I6] ) 

Wliafi = xT^WlT'al (2) 

for that approaches to exp(— 1^) as (Euclidean time) T — >■ 00, and then subtracting the energy of 
two separated qq clusters from the Vi to get the lattice-generated qqqq binding energy Ei for the data 
point. The values of the \^-^\ and thus of Ei, depend only on the numerical values of the elements VF^-^j 
of the matrix of the Wilson loops. The values of k and / depend on the number of Wilson loops evaluated; 
for the qqqq system these were taken to be 1 and 2. Two of the corresponding Wilson operators (whose 
vacuum expectation values are the Wilson loops W'^i) are shown for example in Fig. 1.5 of [M]. Knowing 
the Wilson loops, the procedure of getting Vi can be found for example in eqs. 4, 11, 12 and 15 of ref. P^ . 
The arguments for continuum limits being achieved before extracting Ei are given in ref. |17| . 

The Mi in eq. ([T]) are obtained by subtracting the energy of two separated qq clusters from the eigen- 
values of a matrix obtained through a model of the qqqq system. For this, the model has to give a basis 
and an operator V whose representation with respect to the basis gives a potential matrix V . Mi are 
obtained by setting the determinant of ^ — (Afj -I- Vii)N equal to zero, with Vn being the energy of 
two separated qq clusters and N the (overlap) matrix of an identity operator in the basis. Searching 
for the model, the simplest way to extend a two-particle potential model to a few-body is to use the 
potential for each pair of particles in the few-body system and simply add up such two-body potentials. 
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This approach has been successful in atomic and many-nucleon systems; the corresponding two-body 
interaction being described by Coulombic and Yukawa potential, for example. For a hadron (or a system 
of hadrons) composed of many quarks, antiquarks and the gluonic field, the lowest order perturbative 
Feynman amplitudes are of this sum-of-pair-wise form. Though Feynman diagrams themselves become 
impractical for typical hadronic energies because of larger couplings, models have been tried which simply 
replace the two-body Coulombic potential (essentially a Fourier transform of the Gluonic propagator) by 
more general Coulombic-plus-linear form; see ref. |18j . This approach is not free of problems; for example 
it leads to inverse power van der Waals' potentials [12] between separated colour-singlet hadrons which 
are in contradiction with experimental data. But this model has many phenomenological successes and 
it is worthwhile inquiring if 

1) it provides a basis and operator to generate a potential matrix, and 

2) how good is the chisquare if the eigenvalues of the resulting matrix are used as Mi in eq.([T]). 
The answer provided by ref. and earlier related works is that the model does generate a matrix of 
the required kind. But the resulting chisquare, defined by eq.Q, is too bad; see Fig. 4 of ref. [20]. To 
refine the model we can improve the basis beyond the {|1), |2), |3)} defined as [HI HI] 

|1) (gi93)(g2g4), |2) = ((7194) (9293), and |3) = {qiq2){q3q4), (3) 
and the operator beyond 

4 ^ 

^ = - E K + 1|] + E "»^F-F„ (4) 

i=l ' i<j 

with Vij being is the potential energy of a qq pair with the gluonic field between them in the ground 
state. Or both the basis and the operator can be improved. What ref. [TT] does is to improve directly 
the matrix (representation) after writing down the underlying basis. The authors do this in a number 
of ways. One model, termed model II, uses the same {|1), |2), |3)} basis but multiphes the off-diagonal 
elements of the overlap and potential energy matrices (that is, the representations of the identity and 
potential operators 

i<j 

respectively) by a few-body gluonic field overlap factor / — exp{~bskfS) with bg as the tension of the 
string connecting a quark with an antiquark, S the area of a surface bounded by external four lines 
connecting two quarks and two antiquark and kf = 0.5 approximately; theoretical arguments suggest 
S should be the area of the corresponding minimal surface, though in ref. [11) half of a sum of four 
triangles was used for numerical convenience. At the quark level, this model II was noted to much 
reduce the chisquare of eq.Q. This model has been worked out in [22 [52 [23] till meson-level transition 
amplitudes. The dynamical calculations require a kinetic energy term as well. As this is taken to be 
apart from some technical considerations of hermicity, proportional to the overlap matrix and hence its 
off-diagonal elements are also multiplied by the overlap / factor. Thus / provides one parametrization 
that connects QCD simulations with hard experiments. 

But model II is not the best model of ref. [11 ; the paper continues to then improve the basis by 
including the gluonic excitations as well. That is, it extends the {|1), |2), |3)} basis by including the 
states 

= (gi 93)9 (92(74)9, |2*) (gi'74)g(<Z2'73)g, and |3*) = (9192)9(9394)9- (5) 

Here (9193)9 denotes a state where the gluon field is excited to the lowest state. (The excited states of 
gluonic field can, for example, be seen in the QCD numerical simulations; see ref.[25l[T0] and others). 
When the overlap, potential and kinetic energy matrices are written in this extended basis, their order 
increases to 6 x 6 rather than previous 3x3. If in addition, we introduce many body terms in this extended 
model, new kind of gluonic field overlap factors (/", f^) appear in the off diagonal terms resulting in 



what ref. [TT] terms model III giving the least chisqure in ref. [TT]; (see eq.(15l below); our truncation 
to 4 X 4 matrices is explained before this equation. As mentioned above, the purpose of our present 
paper is to work out this improved model III of ref. llj till the meson-meson scattering amplitudes and 
energy shifts. As this improved model III includes the gluonic excitations, it consider transitions from 
three ground state quark states to the ones having gluonic excitations. And by adding to it the quark 
motion (wave functions) to reach the hadron level, we are now able to study transitions from ground 
state meson-meson systems to hybrid-hybrid systems. 



3 



A worth-mentioning aspect we have studied is the hadron-level imphcations of the differences of 
the gluonic-excitation-including model III and the sum-of-pair-wisc approach. The quark-level work in 
ref. can be interpreted to mean that with suitable extensions and modifications something like a sum- 
of-pairs approach can be a good approximation for a tetraquark system as well. It states " At the shortest 
distances, upto about 0.2fm, perturbation theory is reasonable with the binding being given mainly by 
the |1), |2), and |3) states inter-acting simply through the two-quark potentials with little effect from 
four-quark potentials". However, ref. [11] states, "for large inter-quark distance (greater than 0.5 fermi), 
quark-pair creation can no longer be neglected. However, in the intermediate energy range, from about 
0.2 to 0.5 fm, the four-quark potentials act in such a way as to reduce the effect of the |1), |2), and 
1 3) states so that the binding is dominated by the |1*), |2*), and |3*) states, which now interact among 
themselves again simply through the two-quark potentials with little effect from four-quark potentials." 
This suggests that models involving only two-quark potentials could be justified provided excited gluon 
states (such as |1*), |2*), |3*)) are included on the same footing as the standard states |1), |2),|3). We 
have checked if such features survive at the experimentally meaningful hadron level, by comparing the 
dynamical implications of 

(1) a model extended to the gluonic excitations but otherwise sharing many features with the sum- 
of-pair-wise approach, with 

(2) a model that includes explicit many-body terms but does not include gluonic excitations. 

Thus we report if after including the gluonic excitations a sum-of-two-body potential model can replace 
to some extent many-body potential terms in a tetraquark system even at a hadronic level. Specifically, 
we have calculated in both kind of models meson-meson transition amplitudes Tij from (i) a ground state 
meson- meson clustering to (j) a different ground-state clustering and to a clustering of gluonic-excited 
mesons. Tij are elements of the meson-meson scattering theory T-matrix, can be termed as phase shifts, 
transition potentials or meson-meson coupling, and their absolute squares give meson-meson differential 
cross sections 18J. Moreover, using these transition amplitudes in the second order perturbation theory, 
we study shifts (in both kind of models) in a ground state meson-meson energy due to coupling to a 
different ground-state clustering and to an excited state meson-meson system (i.e. to hybrid loops). 
These energy shifts are what are also termed polarization potentials [26]. 

To reach the hadronic level, we have included the quark motion through quark wave functions. To find 
the quark position dependence in multiquark systems, a number of methods are used such as variational 
method p71 128j . Born-order diagrams [18 , and resonating group method Variational approach is 
used by Weinstein and Isgur to optimize a meson-meson wave function in a quadratic |27| and later 
coulomb plus linear |28j potential and a hyperfine term combined with a sum-of-two-body approach. 
Then they projected the meson-meson state onto free meson wave functions to estimate a relative two 
meson wave function which gives the equivalent meson-meson potential and obtained the meson-meson 
phase shifts. The similar results can be obtained by using Born-order quark exchange diagrams [18 in 
a non-relativistic potential model to describe low energy scattering of qq mesons. 

In the present paper, we have used a formalism of the resonating method as used in |22| . In the 
resonating group method, the dependence on the internal co-ordinates of the system is specified before 
solving the problem to integrate out the degrees of freedom corresponding to the internal coordinates 
of clusters of the system. At a later stage, because of the complexity of the calculations we use a Born 
approximation to specify the dependence on the vectors connecting the centers of masses of our mesonic 
clusters as well. Moreover, we have not included in our basis an explicit diquark-diantiquark state. In the 
weak coupling limit such a state is a linear combination of the meson-meson states and thus cannot be 
included in a basis. Away from the weak coupling this can be included. But inclusion of a third clustering 
states in the basis did not affect the lattice simulations reported in ref. [TT] . Later works [T31 [H] using 
SU(3)c were done with only two clusterings. Thus we have expanded our two-quark two antiquark wave 
function in a basis that is limited to four meson-meson states: ground and excited states for each of the 
two possible qqqq clusterings. 

In this exploratory work, we have taken all the constituent quark masses to be same as that of a 
charm quark and we have taken all the spin overlaps to be 1 without calculating them. Actually we 
have not studied the effects of the quark degree of freedom, not its orbital angular momentum as well; 
we want to study only dynamical effects of including the gluonic excitations in our basis for the same 
quark configurations. In section 2, the potential model in the extended basis in the pure gluonic theory 
is introduced for the static quarks. Basically, in this section we tell where does the model of ref. [TTj 
fits in our full scheme that incorporates the quark motion through a resonating group method formalism 
meaning pre-specifying quark-antiquark wave functions within qq clusters. The coupled integral equations 
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Figure 1: (a) and (b) describes the meson-meson topologies and (c) describes diquark and diantiquark 
topology. 



for the remaining inter-cluster wave function xk^Rk) are written in section 3. In section 4, these 
integral equations are solved to calculate the transition amplitudes and energy shifts we are studying. 
The numerical results for meson-meson system transition amplitudes and energy shifts in both the above 
mentioned kinds of models with concluding remarks are given in section 5. The partial wave analysis 
results of transition from ground state to excited state gluonic field are also reported in section 5. 

2. q^f POTENTIAL MODEL (IN THE EXTENDED BASIS) 

Using the adiabatic approximation, the total state vector of a system containing two quarks, two 
anti-quarks and the gluonic field between them can be written as sum of product of quark position 
dependence function (ri, r2, rg, rj) and the gluonic field state | k)g. (The gluonic state | k)g is defined 
as a state approaching to colour state | fc)c in the limit of quark anti-quark separation approaching to 
zero. Here fc = 1, 2, 3, 1*, 2*, 3*.) The function (^A:(ri, i"2, rg, rj) can be written as 

(f>K{ri,r2,r^,r^) = 0j<:(Rc, Rk, Yif, zk), 

with K = 1,2, 3. Rc is the overall center of mass co-ordinate of the whole system. 
With the notation of Fig. 1, the relative co-ordinates Ri, R2, and R3 are defined as 

Ri = ^(ri -I- rg - r2 - rj), = ri - Tg, Zi = ra - rj 

R2 = ^(ri + rj- ra - r3),y2 = ri - r^,Z2 = ra - rg 

R3 = ^ (ri + r2 - rg - rj), yg = ri - r2, Z3 = rg - rj, 

Ri being the vector joining the centers of mass of the mesonic clusters (1,3) and (2,4); similarly about 
R2 and R3. Now using the resonating group method, the quark position dependence function can be 
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written as a product of function of known dependence on Rc,y/^,Zif and of unknown dependence on R/^. 
i.e. 

(j)K{'Rc,'RK,yK^^K) = V'c(Rc)xK(R'i^)V'K(y/f,ZK)- (6) 

Thus, the two quarks two antiquarks state vector can be written as 

k 

where 

■0fe(yi<-,ZK) = Cfe(y/f)Cfe(z/f), 

Cfc(y_R') s-^d ^fe(z;f) being the normaUzed solutions of the Schrodinger equation for quadratic confining 



potential (written in eq.( 10 )) for a pair of quark-anti-quark within a cluster. We take, for the zero relative 



orbital momentum {() of a quark w.r.t. the antiquark of the cluster, 



1 

^k{zk) = . „,3 exp(--^). 

Here d is the size of meson (detail is written after eq.( |10P ) and m being the constitute quark mass. In 
our case, m is the mass of c-quark and equal to 1.4794 GeV as used in [5U] . 

After writing the form of the wave vector, we describe our Hamiltonian, starting with the limit when 
each gluonic field overlap factor / = /" = /'^ = 1. In this limit, the Hamiltonian whose representation 
matrices in the basis {|1), |2), |3), |1*), |2*), |3*)} would become those in ref. [TT] is 

4 ^ 

1=1 i<J 

where Vij is the potential energy of meson for the ground state gluonic field and Au*j is the difference 
between ground state and excited state gluonic field potential. We take the kinetic energy in the non- 
relativistic limits. This limit is also used in a recent work by Vijande ref. |31j that deals with multiquark 
system (two quarks and two antiquarks) to study the spectrum using a string model for the potential. 
In ref. [11, potential energy matrix elements are written so that the potential energy for each pair ij is 
equal to Vij for the matrix elements of the Hamiltonian between gluonic ground states, and it is equal to 
Vij + Au*- for the matrix elements between the gluonic-excited states. We have modeled these two forms 
by taking e = for the ground state matrix elements and e = 1 for the elements between gluonic-excited 
states. For the elements between ground and excited state gluonic field, the e value that results from 
their parameter ao being fitted to 4 (in their Table 1) is surprisingly 2 and not any value between and 
1. A possibility is that this is a result of them taking the area S, we mentioned above in our introduction, 
in the form of average of the sum of triangle areas instead of the theoretically motivated minimal surface 
area. Thus we have somewhat explored f- — \ between and 1 and e = 1 in addition to e = 2 which we 
have mainly studied. 

In the above equation, Fj (operating on it^ particle) has eight components F'- with 1 — 1,2, 3, , 8. 

Each component is equal to -f, where are the Gell-Mann matrices. I is used in the superscript to avoid 
the possible confusion with index i. We used the potential with the colour structure of one gluon exchange 
in the form given in ref. [11]. With the use of ground state potential Vij in the realistic coulombic plus 
linear form, it becomes impossible for us to solve the integral equations appearing below in eqs.( 21p4 1. 
Therefore we used the parametrization of the static pairwise two quark potential as 

= Crjj + C, with i,j = 1,2, 3, 4. (10) 

In this simple harmonic potential, the parameters C is related to size (d) of wavefunction {^k) through the 
relation C = "i^y = — ^^mo;^ with w = and for consistency of the diagonal term of the integral 
eq.([2l|) C = I (4m + 3w - 2M) GeV [12] with Af = 2m = 2.9588GeV being the mass of a charmonium 



meson. The parameter d is chosen in such a way to reduce the error resulting from a use of this quadratic 
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Figure 2: Overlap of the (ground state) wave function of the reahstic Coulombic plus linear potential 
(shown by solid line) and that of the quadratic potential (shown by dashed line) . 



potential instead of the realistic one. The error may be both in the wave functions of the distance between 
a quark and antiquark within clusters and those of the vectors joining the centers of masses of the two 
clusters. As for the first dependence, we found that the maximum of the overlap integral of each of the 
wave functions S.Kiyx) ^^id (z^f ) of the quadratic potential and that of a more realistic coulombic plus 
linear potential {—^as + bgr) is 99% at c? = 1.16GeV~^. (For the parameters of the realistic potential 
we used values — 0.5461 PO] and bs — 0.1425GeV^ [3D] for mesons composed of charm quarks.) This 
overlap is shown in Fig. 2. A similar work was done in ref. |32j for lighter quarks. They found that the 
overlap wave function of SHO (quadratic potential) and that of coulombic plus linear can be made as 
large as 99.4% with the suitable adjustment of parameters. 

For the additional term in the qq potential for the gluonic excitation, the usual flux tube (vr/r [7]) 
or string based analytical expressions become impractical for us, as mentioned in the introduction. Thus 
for that we tried an ans'dtz of the form of 

Aw* = Ae--^'-'^- . (11) 

This gaussian gluonic potential (^exp(— Brf^)) is a smeared form of as written in appendix of 

ref. [55]. From the Fig. 3 of [TD], we get the potential energy difference between ground and excited 
states for different r values (rj). We choose A and B for which becomes minimum, is defined as 

X^=Y^{e,-AcM-Brl])\ 

i = l 

with n being the number of data points. This gives 

A = 1.8139GeV, B = O.OGSTGeV^. 

For finding the wave function corresponding to our total potential Crfj + C + Ae^^^^^ , we used the 
variational method with an ans'dtz wave function 

CAYk) = "2/i'exp(-ro|.). (12) 
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Figure 3: Overlap of the (excited state) wave function of the realistic Coulombnic plus linear plus 
Aexp{—Br^) potential (represented by solid line) and that of the quadratic plus y4cxp(— i?r^) poten- 
tial (represented by dashed line). 



The normalization of this ^^(y-R') w.r.t yK gives 

n — (422p3)(1557rJ). 

This leaves us with one variational parameter p chosen to minimize the expectation value of the two body 
Hamiltonian in the excited state gluonic field wave function. This gave p — 0.048GeF^. For this value 
of p, the overlap of wave function of the quadratic potential plus Aw*- and that of coulombic plus linear 
plus Ai;*^- within a hybrid cluster became 99.9%. Both wave functions are shown in in Fig. 3. Having 
much reduced the errors in the in-cluster factors of the total wave function, the question remains how 
much the inter-cluster factors of the (terms of the) total state vector are affected by our use of convenient 



but not realistic qq potentials. For the inter-cluster wave functions, eventually we use below in eq. ( 25 ) 
plane wave forms which get their justifications from the validity of Born approximation for our problem 
regardless of potential expressions we use. This plane wave form has only one usual parameter (the wave 
number) and eq.(ClO) below relates its value for the ground as well as excited state inter-cluster waves to 
the very good values of d and p that almost give realistic ground state and excited state wave functions 
within qq clusters. But the relations between the inter-cluster wave numbers and the d and p assume a 
quadratic confinement and this may affect our numerical results but hopefully not at least the qualitative 
features we are pointing out. Perhaps it is worth mentioning here that properties of q^q^ systems were 
calculated using quadratic confinement in ref. [27], and then with the realistic potential in ref. |28j and 
both the works favoured the existence of meson-meson molecules. 

Now we combine our Hamiltonian and all the wave functions we have mentioned in the Schrodinger 
equation for the meson-meson system, which means that the overlap of (H — E^) \ ^) with an arbitrary 
variation | ^^E*) of state vector | XP) vanishes where | is the state vector of the whole q'^lf system. In 
I 5"^), we considered only the variation in Xk (see eq.(l7])), as in the resonating group method. Thus we 
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wrote 



(5* I H - I *) = ^ / d^R,d^RKd''yKd''zKM1<^c)6xki'RK)UyK)U^K) 

k,l 



(13) 



forfc,/ = 1, 2, 3, 1*, 2*, 3* and K,L = 1,2,3. The arbitrary variations 5xfc(R-K)'s for different values of 
Rk are linearly independent and hence their co-efficient in eq.(13l should be zero. With the trivial Rc 
integration performed to give a finite result, this leads to 



^ J d'yKd'zKUyK)U^K)g{k \H-E,\ l)gXl{^LMyLM^L) = 0, 



(14) 



where 



j{k\Il-E,\ l)g ^g{k\KE + V + Am-E,\ I), 



Elements of V and KE matrices are defined below in eq.(17l and eq.(18|. In ref. jTT] it is stated that, 
in the lattice QCD simulations, it was found that the energy of the lowest state was always the same in 
both a 2 X 2 and 3x3 description, provided | 1) or | 2) had the lowest energy. In addition the energy 
of the second state was, in most cases, more or less the same. Two level approximation is also used in 
later works [131 [H] of the tetraquark system. Considering this, we include only two topologies (1 and 2), 
meaning four states(| 1), | 2), | 1*), | 2*)). According to the model of ref. [IT] the overlap matrix in this 
truncated 4-basis is given by 



N = {g{k\l)g}- 
And the potential matrix V{f) is 

V^{g{k\V\l)g} 

Here, 
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^Ks + ^24) 



V12 = 1^21 g/(t'l2 ^13 - 



- '"li - ^^23) 



^^22 = -3(Vl4+«23) 
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V21* = = - 

Vi.2 



V12* = - 
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2(<3 + "24) - ^(^J'i + ^^23) - ^(-«?2 - ^'li)) 



Vi.2* = V-2n* - - Y^/M - «3 + ^^24 + <4 + ^"23) + 10(^'^2 + V^) 



V2*2* 



1 

' 6 



6K4+"23) 



V2- 



(17) 
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with v*j = Vij +eAv*j, e being defined above (after eq.(|9|)). The coefhcients of Vij and w* , resulting from 
the F.F operator, are given in Table 1 in the Appendix A. The kinetic energy matrix of the two quarks 
two anti-quarks is taken to be 



4 

KE = {,{k I KE I I),} = N(f)l{—Y,Vl)N(f)l. 

i=l 



(18) 



The kinetic energy in the same form is also used in ref. [22]. 

The gluonic field overlap factor /, as written in introduction, is suggested by ref. ^11^ as 

f = eM-b.kfS], (19) 

with bg = O.lSGeV^ [7]. In ref. [13] the gluonic field overlap factor / is used in the Gaussian form as 

f = exp[-kfbsJ2rU (20) 

i<j 

employed in SU{5)c for interpreting the results in terms of the potential for the corresponding single 
heavy-light meson. In ref. |I3j, the simulations that are fitted by using / are for the configurations when 
the gluonic field is in the ground state i.e. overlap matrix is a 2 x 2 matrix. In ref. [11) . simulations 
are reported with 2-colour approximation. But in ref. |13| . lattice simulations performed for SU{^)c are 
reported. Our overlap, potential and kinetic energy matrices are also written in SU(3)c, so we use the 



kf multiplying sum of area form of / (written in eq.(20)) with fc/ = 0.6 (as used in |13j ) for numerical 
convenience and have not used the minimal area form. When we are observing the dynamical effects 
for the ground state, our overlap, potential and kinetic energy matrices are 2x2 matrices and we use / 
with kf ^ 0.6. But when we incorporate the excited state gluonic field our overlap, potential and kinetic 
energy matrices become 4x4 matrices, in the upper left 2x2 block of these matrices, the form of / 
remains the same but the value of fc/ is changed to 1.51 according to the conclusion of ref. [IT]. In the 
other blocks are also used and defined in refs. [IT] [3T] as 

r^{f{' + bJSS)exp[-b,kaS], 

f = exp[-6sfccS']. 

If we take /" as a function of area as defined in ^llj, it becomes unmanageable to solve the integral 
equations ( 2l]|24 | and hence we have taken /" to be a constant but have tried a variety of its values to 



explore how much our conclusions depend on its value. As for f^, the fit in ref. [3T] [TT] of the model to 
the lattice data favours fcc = which implies that — 1 i.e. the excited configurations interact amongst 
themselves in the way expected from perturbation theory. Thus we have used f = I. 

3. COUPLED INTEGRAL EQUATIONS 



Using "A^"," 1^" and "A'i?' elements in eq.p^, we got four integral equations for four different values 
of k or I. Then we do and zx integrations. All the integrations required are in the Gaussian form 
or modified Gaussian form (with a polynomial in the integrand multiplying the Gaussian exponential) 



and we integrate analytically. For K = L = 1,2, in eq.(14|, x/(Ri) is independent of and zk and, 
thus, can be taken out of integrations. After the integration, the result is R/f dependent co-efficient of 
Xki^x)- For K ^ L, and are replaced by their linear combinations with one of them as identical 
to Rl and other one independent of it as R3. The jacobian of transformation is equal to 8. Then we 
integrate the equation w.r.t R3. Integration leaves the following four equations: 

Ai(Ri)xi(Ri) + y d3R2i^(Ri,R2)x2(R2) + y"d3R2£;i(Ri,R2)x2(R2) = 0, (21) 
A2(R2)X2(R2) + j d3RiF(Ri,R2)xi(Ri) + j d'Ri£;i(R2, Ri)xl(Ri) - (22) 
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F4(Ri)xt(Rl) + J d3R2i;2(Rl,R2)X2(R2)+ J d3R2i;3(Rl,R2)X2(R2) -0 
f4(R2)X2(R2) + J d3RiS2(R2,Ri)xi(Rl)+ J d3RlS3(R2,Rl)X^(Rl) =0. 



(23) 



(24) 



The symbols are defined in the appendix B. We have eventually replaced ri, r2, r3, r4 by Ri, R2, R3, and 
Re- With trivial integration on Rc, we have eq.( 14) that is independent of Rc 



Now, after the integration 

on R3, the above four integral equations ( 2Tp4 | depend only on Ri and R2. So every quantity which we 
want to calculate depends on Ri and R2. In eq.( 2l]p2 1, the first two terms in each equation, containing 
Xi(Ri) and X2(R2), are for the ground state. It is noted that in these terms (observing the definitions 
of the symbols in appendix B), there is no dot product of vectors Ri and R2. So the results from these 
terms should not depend on the angle between Ri and R2, called 9. The third term in each of eq.( 2Tp2 ) 
is due to the gluonic ground and excited states. In these terms dot product of two vectors (Ri and R2) 
appear, so the results from these terms depend on 9. 



4. SOLVING THE INTEGRAL EQUATIONS 



Now taking the three dimensional Fourier transform of eq.( 21|23 ) with respect to Ri and eq.( 22|24 ) 



with respect to R2, we get formal solutions xi(Pi), X2(P2), Xi(Pis), and XiCPis) as shown in appendix 
C. Because of the coupling to the gluonic excitations, it become difficult to solve the integral equations 
for non trivial solutions for Xi(Pi)j X2(P2), Xi(Pis)' X2(P2s) analytically as done in [551 [53]. In 
[55] , the meson wave functions, including the gluonic field overlap factor, is Ri, R2 separable. So there 
the integral equations can be solved analytically by replacing xi and X2- But in our present work, the 
meson- meson wave functions are not separable in Ri,R2. So we use the Born approximation (as used 
in [18 for meson-meson scattering) to solve the integral equations. Our results given below also justify 
our use of the Born Approximation. Using this approximation, we use the solutions (xi(Ri), Xi'(Ri)) of 
eqs.( 2Tp4 ) in absence of interactions (meaning / = = /'^ = 0) 



Xi(Rz) = \ - exp(zP,.Ri), 

[2 ^''^ 
and x*(Ri) = Y- exp(zPi5.Ri) 

for 1 = 1,2. Here the coefficient of exp(iPi.Ri) is chosen so that it makes x(Ri) as Fourier transform of 
^^^^2^^y^. Similarly the coefficient of exp(zPis.Rj) is chosen. Using this approximation, the integration 
on Ri and R2 can be performed to get xi(Pi) (written in eq.(Cll)). 

Til can be calculated (as in ref. j22j) by considering the coefficient of /^^^p^^ containing the xi(Ri) 
from eq.(Cll). As in this equation, there is no coefficient having xi(Ri), so it gives Tn — 0. T12 can be 
calculated by considering the coefficient of ^^^p^^ containing the X2(R2) from eq.(Cll) in the following 
eq. 

T12 = Af JPi[coef. of ] containing X2(R2)], (26) 

with M being the mass of cc meson. Similarly T12* can be calculated by considering the coefficient of 
^^^P^^ containing the X2(R2) from eq.(Cll) 

ri2* = A/^Pi,[coef. of ] containing X2(R2)]- (27) 

The relation between off-diagonal transition and scattering matrix element is written as 

S = 1 - 2lT, (28) 



where 5', T, and I represent 4x4 scattering, transition, and identity matrices respectively. The eq.( |28[ ) 
can also be written as 

Sij = 5ij - 2lT.ij. (29) 
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for i,j = 1, 2, 1*, 2*. Using the transition matrix elements, the contribution to the energy shift of meson- 
meson system (cccc) through £ — states (both with and without gluonic excitations) can be calculated 
by using the stationary state perturbation theory, i.e. 

= + T„; + J2r Jf^dP,, (30) 
Jo E. ~ 

with the initial state i and intermediate state j. We have considered initial states where the gluonic 
field should is in ground state, so i = 1,2, but intermediate j — 1,2, 1*,2*. Here Ta = {i\T\i) , — 
|(j|T|i)p, is the energy of a ground state (1 or 2) of meson-meson system, and may be the energy 
of the other meson-meson ground state or that of an gluonic-excited meson-meson state. 



5. RESULTS AND CONCLUSIONS 

l-The transition amplitude T12, from one meson- meson ground state to other, is calculated by using 



eq.(26) with fc/ =0.6 [TH] (without the incorporation of gluonic excited states). Its dependence on the 
center of mass kinetic energy is shown in Fig. 4 below. As it is noted that magnitude of transition element 
T12 is less than 1, so this result shows the validity of Born approximation. In result 3 we compare the 
transition amplitude of this many body ground state gluonic field model at kf — 0.6 with the transition 
amplitude obtained from a model that is extended to gluonic excitations along with changing kj = 1.51. 

2- The transition matrix element T12* , for transition from ground state to excited state gluonic field 
with kf = 1.51, depend on the parameters e, and /°. T12* also depends on 9 (the angle that Pi makes 
with P2 and Pis makes with P2s)- We take parameter as a constant as discussed earlier in section 2. 
For e = 2, we take different values of to see the effects of /" on T12*. Fig. 5 shows the dependence 
of T12* on f"" at 9 = 90. By taking different values of 9, T12* was calculated. We found that the results 
are slightly different for different 9. This slight angle dependence is not directly reported here, but can 
be found by a linear combination of the corresponding m = spherical harmonics with coefficients for 
each value of energy read from Figs. (8-10) that report the partial wave amplitudes that result from 
this angle dependence. For e = 1/2, the dependence of on is shown in Fig. 6. And for e = 1, 
the dependence of T^2 is shown in Fig. 7. These graphs show that the magnitude of transition 
amplitude is increasing with the increase of 

3- For the center of mass kinetic energy range to 2 GeV, the average of the modulus of transition 
amplitudes (excluding excited states) is equal to 0.0369 at kf = 0.6. But if we change the value oi kf 
from 0.6 to 1.51 and also include the excited state gluonic field, the average meson-meson transition 
amplitude again becomes (0.0218 + 0.0152) = 0.0370 by choosing f = 0.03 and 6* = 90 with e = 2. This 
means that a sum-of-two-body potential model can replace to some extent many-body potential terms 
in a tetraquark system by including the gluonic excitations. (An exact sum of two body terms would 
required /a = 1 though). In above, 0.0218 is the transition amplitude due to ground state gluonic field 
at kf = 1.51, and 0.0152 is the transition amplitude due to gluonic excited states. 

We have also explored the results with e = 1/2 and e — I. As e is proportional to the potential matrix 
elements taken between ground and gluonic excited states, when we replace e = 2 with e = 1/2 or 1 , 
it has no effect on ri2 . But for e = 5 with inclusion of gluonic excitations, now the values /° — 0.08 and 
9 — 75 almost restore the average transition amplitudes for kf — 1.51 to (0.0218-1-0.0149) — 0.0367 . (The 
significance of the first and second term is as before). For e = 1, at = 0.055 and 9 = 90, the average 
transition amplitudes for kf = 1.51 becomes (0.0218-1- 0.0150) — 0.0368. This almost restoration again 
indicates that perhaps we can always choose parameters etc. so that with an inclusion of the gluonic 
excitations the sum-of-two-body potential model can replace many-body potential terms in a tetraquark 
system. 



4- The effects of interaction are also observed on a contribution to the energy shift by using eq.(30). 
Here Ta = with i — 1,2. In J._^^ e9-e° have chosen E^ — 5.9176 GeV, the value of energy for 

which center of mass kinetic energy becomes zero. The intermediate energy state E^ we study depends 
on the state of gluonic field ( |2), 12"^) are only possible states that couple to |1)). For the gluonic ground 
state (|2)), E^ = 5.9176 + 0.3380^^ (calculated using eq.(ClO) with P^{1) = P}), and for the excited 



12 



0.1 




Tc (GeV) 



Figure 4: The graph for the energy vs. T12 for ground state at kf = 0.6. 
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Figure 5: The graph for the energy vs. T12* at ^ = 90 with e = 2 for different values of The curve 
with dots is for /° = 0.015, with lines plus points is for /" = 0.05, with lines is for /" = 0.08, with points 
is for f°- = 0.1, and with steps is for = 0.15. 
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Figure 6: The same as Fig. 5 but with e = 1/2. 
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Figure 7: The same as Fig. 5 but with e = 1. 
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state gluonic field (|2^)), E° = 5.4638 + 0.3380Pf^ (calculated using eq.(ClO) with P^ils) = PfJ. At 
kf = 0.6 and Ec = 5.9176 (i.e. at the threshold), the shift to the ground state meson-meson energy is 
found to be Ei — Ef — 0.7268 GeV excluding the gluonic excitations. With increasing kf to 1.51 as usual 
and including the coupling to the gluonic-excited meson-meson state, we want to restore the same energy 
shifts we can get to (with same origin of the first and second correction terms) 



E, = K° + (-0.0950 - 0.6371)GeV = E° - 0.7321GeV 

if /° = 0.14 is used. This shows that gluonic excitations can replace the many terms for £ = energy 
shifts as well. We note that the energy shift is independent of 6, the angle between Ri and R2. 

The energy shifts we have reported here can be compared with the hadron-loops-generated mass shifts 
to charmonium states reported in Table III of ref. |33| . It is difficult, though, to conclude anything from 
this comparison as in ref. [53] the integrand contains squares of the matrix elements of the ^Pq meson 



decay amplitudes whereas in our integrands in eq.(30l contain squares of meson-meson couplings. Thus 



though the intermediate states in both works are respective hadron loops qqqq, the initial and final states 
in ref. [33] are qq but in our work initial and final states are also qqqq. Only we include the glounic-exited 
intermediate qqqq states (i.e. the hybrid hadronic loops) for our problem. 

In result 3, it is noted that the average transition amplitude obtained (for e = 2) by a model that does 
not include the gluonic excitations is equal to the average transition amplitude obtained by including 
gluonic excitations for /° = 0.03, but the energy shift obtained by both models becomes comparable at 
/° = 0.14. One possible reason of this difi^erence in the values of /" could be that the average transition 
amplitude is calculated for the center of mass kinetic energy range in between GeV to 2 GeV, but 
energy shift is calculated at threshold center of mass kinetic energy. 

5- For the ground state gluonic field, transition matrix elements from state 1 to 2 and from 2 to 1 are 
9 independent. But for transition elements to the gluonic-excited meson- meson state depends on 9. We 
projected this angle dependence on spherical harmonics Yim{d)- The results of this partial wave analysis 
are also reported in Figs. (8-10) for m = and /° = 0.03 (this value is used above in result 3). The 
reason for truncating the spherical expansion to to = harmonics is that we have no dependence on 
4>, the azimuthal. These analysis shows that partial wave amplitudes are decreasing as we go from Yoo 
the coefficient to the Yso coefficient i.e. from S-wave to H-wave. Figs. 8 and 9 are for the even wave 
{D,G,...) ratios with S'— wave. S/D ratios are also used in ref. [34l[32]. These graphs shows that the 
S/G ratio is too much large as compared to the S/D ratio. It means that 5— wave is dominant over 
G— wave. In Fig. 10, ratios of odd waves with 5— wave are shown. We noted that S/H is too much 
large as compared to S/P. This shows that H,J,... waves can be neglected as compared to S"— wave. 
The partial wave analysis indicates the presence of P, Z?, P, G, H waves only when we include the gluonic 
excitations in combination with essentially sum-of-pair-wise approach. It means that, in the presence of 
gluonic excitations, an = 0, 1, 2, 3, ... ground state meson-meson system may couple to £ — 0, 1, 2, 3, ... 
hybrid-hybrid systems as a intermediate states or as final states. 

As mentioned above in result 3,4, we compare our meson-meson transition amplitudes and polariza- 
tion potentials (in the respective center of mass energy ranges) in an extended almost sum-of-pair-wise 
approach and a many-body-term based one. Apparently this comparison has even a qualitative differ- 
ence in case of the transition amplitudes {T12 and T12) which have a dependance on the scattering angle 
(meaning coupling to P and higher orbital excitations in the final channel) if we use gluonic excitation in 
combination with essentially sum-of-pair-wise approach (can see Fig. 5) but we get the same transition 
amplitude for all the scattering angles if we use a many-body potential linking ground-state clusters only. 
But a recent work [24J indicates that the absence of angle dependence is an artifact of the unjustified 



overly simple form of the gluonic overlap factor / (defined in eq.(20l); using a proper lattice-gauge theory- 



based form of / (defined in eq.(19)) also results in this angle dependence and the resulting coupling to 
higher orbital excitations. 

The meson-meson to hybrid-hybrid coupling worked out in this paper may affect properties of any 
system that is presently understood to be a purely ground-state meson-meson system. This highlights, in 
addition the energy shift of the meson-meson molecules due to coupling to hybrid-hybrid systems, pointing 
out the possibility of a hybrid-hybrid component in the wave functions of mesons like X{3872), ^(3940), 
y(4260), and Z(4433) which are considered to have q^lf components. And in the above paragraph, we 
indicate that this coupling may also include coupling to orbital excitations. Thus we hope to point out 
a new dimension in discussions about the structure of some exotic mesons as well. 
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Figure 8: vs. S/D ratio. 
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Figure 10: vs. S/P (Dashed line), S/F (thin sohd line), S/H (thick sohd line) ratios. S/P and S/F 
are so small as compared to S/H that the curves for S/P and S/F are lie along horizontal axis. 
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Table 1: The matrix elements of the Fi.Fj operators 





< 11 11 > 


< 11 12 > 


< 2 .. 2 > 


< 11 12* > 


< 21 11* > 


< 1*1 12* > 


< 1*1 11* > 


< 2*1 12* 


FI.F2 
FI.F3 

F1.F4 

F2-F3 

F2.F4 
F3.F4 
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4 
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10 
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APPENDIX A 

Table 1 is used for the matrix elements of the Fi.Fj operators (introduced above in eq.(|9])), for 
different values of indices i and j. In this table ground state elements (with out star) are taken from [22] 
and excited state elements (with star) are calculated by using the following group theory results: 

||3l2334 > +y'^|6i2634 > 



II13I24 \/^|3i2334 > +^^1612634 > 



I814823 >- \l 9II13I24 > ~\l g|8l3824 > 



II14I23 >= \/^|ll3l24 > +a/ ^1813824 > 



3, -..-.4 ' ■ V 3 

and 



|3l2334 > +\ ^1612634 > 



< 8i3824|Fj.Fj|8i3824 >- -. 



APPENDIX B 

The terms used in eqs.( 21p4 ) are defined as 



6 Zm 2 



rm ^ y / -(l + |fc,M^)(Rf + R^) \ i;, 8Cd 

^ Urf'(l + ifc/M2); V )^ 3 3(l + |fc/6,d2) 



I*- + 7r--H ^ - (R-i + «'2)( 



3 9 6md^ V 2 ^ 1 ' 
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n2,R? + 80/Rf ) + __^(400.^R? - 840p) + ,(|^}, 
i;(Ri) = exp(-fiR2)exp(^4?V^)' 



i^2(Rj,Rj) = 2R^ 



+ 



2d2 

10B2r2 2B4Rt 15 4B2r2r2 

J I J I I '■J 



2(R,2) 8S2(Ri.Rj)2 



£;(Ri,R,) = {^-^^-—)Hl6EiKj)F2{Ki,Kj) + 2EiKi)F2{Ri,Ri)}, 
2^ + 2p + i? 



+ i^"''? - + - + ^ + SCST^'"' " ^ ''''''' ~ ^ 

+ ^ + + 4(37^2^)2 (-^2^+ ^^^^^^ + ^ + V) + ^(.^ 



d2 #^2(2i.+2p) #(_i^+2p)2J' 

4r225Sc. TT 3 _ ttV 225A 225c 1575c \ 8 / 105 31_ 

i^4(Rij-4m + n l-^(-) 

15 299 3 685 39690 2 10395 315 6615 \ tt 3 8V^^ 2257r3 



64p2i6p 8p256p2 8p3 (8p)5 p3 (8p)V Mp' 2m 16384p' 

F5(R., R,) = R^R| + + + 15R| - 36R2r2 - 32(R,.R,f) + ^(2R|R? 

F6(Ri,R,) = (-10Aexp(-BRi + r/) - 10^exp(-SRi - r/)){^(r4 + + 4R2r2) 

T54t>4 6 6 t>4t52x 210 2 „2x 945 R^R? WOR^ R^R^ 
+ ^i^J + ^i^i + ^i^j) + (^(^' + ^i) + (8^ - - "(8^ - ^^"(8^ 

_ R^ _ 140R| 32(R,.Rj)^ 16fl|i?| 
2p (8p)3 (8p)2 ^ (8p)2 ^' 
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i;i(R,,R,) = ^^^^^^^^^ 3 cxp( ^ )exp(-2pR,^) [( ) ^ ^^i (R» , Rj) - 2AE(Ri,Rj)], 
E2{B.i,Rj) = ^^-^^^ J d3R.exp(^)exp(-2pR2)[(_^^^)ii^3(R„R,.)-2A£;(R2,-,R,)], 



E3{Iii,-Rj) =— ^exp(-2pR^ - 2pR;,. 



( J ) i (i;,F5 (Ri , R,- ) + Fe (R^ , R, ) ) + 2^ 



3 

[exp(-BR2){^(R„ R,)i?i(R,) + RfR| + 2(R,6 + IitR^^)D2{n,) + 2(r2 + r2)Z?3(Rj)+ 
L'4(Rj) - 4R|i?5(Ri, Rj) - 41?6(Ri, Rj) - 8R^'L'7(Ri, Rj) - 4R^I?8(R,) - 4I?9(R,) 

- 8R^Dio(R,) + 16£)n(Ri,R,)}exp(^^|^) + exp(-BR?){£)(R,,R,.)£)i(Ri) + R^Rj 

+ 2(R^ + R^R^')£>2(Ri) + 2(R? + R|)r'3(Ri) + £>4(Ri) - 4R^£>5(Ri, Ri) - 4£'6(Ri, Ri) 

- 8R^D7(Rj, Ri) - 4R^£)8(Ri) - 4£>9(Ri) - 8R|£>io(Ri) + 16£'ii(Rj, R^)} 

, (2B)2r2 1 , tt ,3 _ 

^^p(i(^)] + 2;^(4^)'^-(^-^^ 



L»(Ri, Rj) = Rj + R- + 4R2r2, 



10(25Rj)2 (2SRj)4 15 
~ {2B + 839)3 (2B + 833)4 (25 + 8p)2 ' 



(2BR,)2 



{2B + 8p)2 {2B + 8p) ' 

(2.BR,)6 21(2BR^)4 105 105(2ER,)2 

(25 + 8p)6 + (2B + 8p)5 (25 + Sp)^ ^ (25 + 83))* ' 

(25Ri)8 36(25Ri)^ 378(25R,)4 126G(25R,)2 945 



£)4(Ri) 



(25 + 833)8 (25 + 8p)7 (25 + 8p)6 (25 + 83>)5 (25 + 833)4' 



(25 + 8i>)2 (25 + 83))' 

(25)6(R,.Rj)2r4 ^25)4 2 2 2 4 (25)^ 2 

^^(^^'^^■^ = (25 + 8p)6 + (25 + 8p)5^^^^-'^^*-^^'^ + ^^^^^ + (25 + 8p)4(^^^^*-^^"^ 
+ 14R?R| + a^a^ + + + ^^j^, 

2 



- JWTW^''^^"^'^' ^ (25 + 8p)3^^(^*-^^^' + ^^^^'^ + (25^?8p)^' 



D8(R,) = 



_ (25)2R| R| 
~~ (25 + 8p)2 25 + 8p' 
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^ 35R,^ i2Br s I 19(2i?)^ e , 77(2i3)^ . 

{2B + 8p)3 (25 + 8p)6 (2B + 8p)5 ' {2B + 8p)4 J' 



(25)4 (25)2 5 2 

(25 + 8p)4 ^' ^ (25 + 8pf ^ ^ (25 + 8p)^ 



5»ii(R„Rj) - ^^^-p^-yjR|(R,.Rj)^ + ^^^--p^^(5Rj(R,.Rj)^ + R-R|) 
1 



(25 + 8p)- 



,(2(R,.R,)^ + R,^R,^), 



5'i2(R^, Rj) = ^^^^ I - 9975 + 32768/R^Rj*(5,2 + R^^) + 40p(287i?,2 + 3ii?2) 
+ 4096p^(2i?fi?| - 43i?fi?| + 2i?2i?6) 

+ 512p3(i?2 + Rl)(l5Rf - 14i?2i?2 + 15^4 ^ i6(2(R,.Rj)2 + r2r2)). 
64p2 ^2875^^ - 772RjR] + 287i?| + 672(2(R,.Rj)2 + R- R^) 



- 1024(7R,R| + 14(R,.Rj)' 

APPENDIX C 

Solving the Integral Equation 

By taking the three dimensional Fourier transform of eqs.( 21|2"3 ) with respect to Ri and eqs.( 22|2"4 ) 
with respect to R2, these integral equations become 



(-5, - 8C<f -\c+^ + \uj + 4m)xi(Pi) + f-J^TTvVrT^ ) / rf'Rirf'R2exp(zPi.Ri 

-(l + |fc/b,rf^)(R^ + R^) \ Ec 8C(f ^_^(^ 1 + /15 



,2 ./X2(R2) + ^ /cZ3R,d3R^e,p(,p^.R^)e,p(_|2)^^p(_2pRL^^^^^^^^^3 



[( , )^fi(Ri,R2)-2A5(Ri,R2)]x^(R2)=0, 
232 +2p 

(CI) 



(-i^e - 8Cd^ - 3^ + 3 + 2- + 4„.)x.(P.) + ^ U^-(i + i^/M-) ,^ / d^n^d^n^eM^^.n.) 

-(l + |fc/ b,d^)(R^ + R^) \ 8C/ 1 + /15 

2d^ / " 3(1 + f fc/6,(i2) 1" " 9^ \Y ' ^"'1 ^ '"2^ 



^ 2.2" ^J [-f- 3(i + ife,M2) + ^-9^ + ^^^y-(Rf + ^: 

(^^^^^))]Xi(Ri) + -^/ d3Rid3R2exp(zP2.R2)exp(:^)exp(-2pR^)^^-|g 



TT 



[( )sfi(R2,Ri)-2Aii;(R2,Ri)]xURi) = 0, 

23^ +2P 

(C2) 
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^ 4r 225£;c, TT ttV 225A 225c 1575c \ 8 , 105 31 15 299 

4m + n i -( — r 1 1 H ( 

, ^(4p)4^2p^ 3 Vl6(2p + B)i(2p)i 16(2p)7 32(2p)V 2m^(8p)3 2 64p2 I6p 

3 685 39690 210395 315 6615 tt 3 SP^ 22577^ \ 



8p256p2 8p3 ^ (8p)5 p3 p4 '^4p' 2m 16384p 

(27r)3 18V2(27rd2)i 



(2^)t 



d3Rici3R2exp(iPi,.Ri)exp(^)exp(-2pR^)[ ]x2(R2) 

J d3Rid3R2exp(iPii.Ri)cxp(-2pR?)exp(-2pR^)[ ]X2(R'2) = 0, 



(C3) 



,, 225Ec,TT,o ttV 225A 225c 1575c \ 8 , 105 31 15 299 

4^77 + n^i -( — y 1 1 -4 ( 

^ (4p)4 ^2p' 3 \l6{2p + B)i{2p)i 16(2p)7 ^ 32(2p)8; ^ 2mM8p)3 2 64p2 I6p 

3 685 39690 210395 315 6615 tt 3 8P|^ 22577^ \ 

8^25fp ~ "(8pF ^ ^^^4^^ ~^ 1638VV^ ^ 

+ TO Xf.n^^ I rf'Ri'^'R2exp(iP2s.R2)exp(-^)exp(-2pR?)[ ]xi(Ri) 

(27r)2 18v2(27rd2)2 J 2d'' 

1 8fn^ J ^3p^^^3j^^^^p(^p^^_j^^)^^p(_2pR2 _ 2pRl)[ ]x^(Ri) = 0, 



{27r)i 



(C4) 



where x(P) is the fourier transform of x(R-)- PijP2 are conjugate to R2 for ground state gluonic field 
and Pis, P2S are conjugate to i?2 for excited state gluonic field. The off-diagonal terms of eq.(C3) 
and (C4) are too lengthy and not used for results, so dots are used to reduce the length of equations. 
The above eqs.(Cl-C4) have formal solutions as [22] 

Xi(Pr) = SiP, - P.mPUD - ^(^( ,,2(,/2,^,^,2) ) ' / exp(.P,RO 

^ ~(l + |fc/&,d^)(R^ + R^) ^^ 8Cct , 4m 8 1+|M^/15 2 



2d2 3 3{l + lkfbsd^) 3 9 6md2 ^2 

1 + _ ^^1^^^ ^^^^3 1 d3Rid3R2exp(zPi.Ri) exp(^) exp(-2pi??) 



^ d2 

8fn2 



-[( , ^ ^ )3Fi(Ri,R2)-2Ag(Ri,R2)]x^(R2), 



18\/2(27rd2)f ''^ + 2p' 

(C5) 

X2(P2) - <^(P2 - Pc{2))/P!{2) - ]p..\s ( ,2n^k ^^2, ) ' / '^'Rirf'R2 exp(zP2.R2) 

^2(^2) (27r)2 V7rd2(l + ^kf-bgd^) J J 



'{1 + pfbsd^){Rl + nl)\^ Ec 8Cd ,4m 8^, l + jkfb^d^flb ^^2 

6^ 

,1 + ifcf6sc^^^, . 1 1 „ . ,-Rl 



2d2 3 3(l + |fc/53d2) 3 9 6md2 

(^^^^))]xi(Ri)-^^-^/d'Ri'^'R2exp(zP2.R2) exp(^) exp(-2pi?2) 



A2(P2) (2^)i 

^^""^ -[(^^)iFi(R2,Ri)-2Ai?(R2,Ri)]xl(Ri), 



18y2{27rd2)f ''^ + 2p' 



(C6) 



1 1 Rrn"' f 
xl(Pi.) = - P,{1s))/P!{Is) - ■ , ^3 ^,,3 / d^Rid^R2exp(zPi,.Ri) 

'^il-r^isj (27r)2 18v2(27ra'^)2 J 

exp(^)exp(-2pi?^)[ ]^^(r^) _ . /^"'./.a / d^Rid^R2exp(zPi,.Ri) ^^^^ 

2a^ Ai(Fis) 3 (27r)2 J 

exp(-2pi?2)exp(-2pi?2)[ 



22 



X5(P2«) = 6{P2s - Pc(2s))/P2(2s) 



1 



1 



exp(^)exp(-2piif)[ ]xi(Ri) - 

exp(-2pi?i)exp(-2pi?f)[ ]xt(Ri) 



A2(P2s) (27r)i 18y2(27rd2)f 
1 8fn^ 



MPhs) 3 (27r) 



y c^^RlC^^R2exp(^P2s.R2) 
7r)2 y 



id^R2exp(iP2s.R2) 



(C8) 



where 



Ai(Pi) = i-Ec - iCd? - |c + + + 4m) - te 



2m 2 



A2(Pi«) =4m + n^{— ^(-)-^ - 



225A 



225c 



1575c 



(4p)4 ^2p^ 3 Vl6(2p + B)i(2p)i 16(2p)7 32(2p)V mV(8p)3 2 



4 / 105 31 



15 299 3 685 39690 



+ 



+ 32p- 



10395 315 6615" 



A TT 3 8(P..)2 22577^ 

64p2i6j3 ' 8p256p2 8p3 ' (8p)5 ' (8^)3^^^^^ 2m 16384p7J 



+ 



(C9) 



and 



Pe{2) = P^{l) = \2m{E^-Am 



m<P 3 



P.(2.) = P.(ls) = ^2m(^. - 4m) - (0.6667C + + ^q^^,,, ^\^y/, )m - 64p. 

Using the Born approximation, the integration on Ri and R2 in eq.(C5) can be performed to give 



(CIO) 



Xi(Pi) = .5(Pi-P,(l))/P,2(l) 

2(1 + Ikfbsd?) 



1 



exp 



Ai(Pi) 
0.1 



1,2 1 1 2/ 1 y / 2-K(f 



1 12l 1 ,3 2/l\2 



exp( 



-Pl<f 



2p + 



1<p- 



4A 



6V2 2p + 



n 3 —pi TT 3 TT 3 

l-)^exp(— ^)(;^;^)2(^^)^ 



2d2 



'8p + 4B''2p + S' 'P2(p)' 



exp( 



PiB_\2 



f.p, _4, 

4P2(P) 



4A 



6\/2'2p + 



— )5exp( 



-P2 



)H- 



TT . 3 



4P2(P) 



2d2 



8p + 4B''2j9 + B' 'P2(p) 



) " exp( 



16^ , TT , 3 , —P? , , TT , 3 ^ 

H ^( — r^)'exp( -)( — )^n4 

3^2 P3(p) 8p 



5+^ 



B2 



-)2exp( 



_p2 

Is 



F3(p) 



4(5+23 





2d2 F3(p) 



(Cll) 



0^3 



Here, 
Oi = — 



8Cd 



CI Im I ^^^^ + ^^^'^-^'^ ( ^^ + - Pid'' 



4md2 



2m# 



(1 + Ikfbsd'Y 1 + |fc/6«d2 
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/ -5Pf 15 15 Pi ^ 

[a{2p + B)3 + 4(2p + 5)2 + 4(2^^ + 2py 4{2p + + 2p) 2{2p + B){^ + 2p) 

( 4^2 X nPu-^_ 3 \ 

1, (2p + B)3 + 4(2p + 5)2 ( + 2p) 1, 4J'| (p) ^ 2F2 (p) 

/ -mS 2lB \ / -2iPfB p p \ 1 

^16(2p + B)n 1 ^ ° V 4F|(p) ^ 16^^4(p) +4F|(p)y' 
- 32.BHP,^ + A, + P,.) (((^^ - + + - 64(^^1^ 



8p + 4B^^ ■ 8p + 4B' ' ■ ''Jm{pf F2{p) 



5P2 15 15 _ P2 



4(2p + P)3 4(2p + B)2 4(2^^ + 2p)2 4(2p + 5)2(2^1 +2p 



2(2p + B)(2i^+2p); V(2p + S)3 4(2p + B)2(2i^+2p); V 4F2(p) ^F^ip) 
I I \ / -2.p2j3 \ 1 , 1 [p, 

^V4(2p + P)3^4(2p + P)2(^+2p)A4(2p + B)'^ 'V 2F2(p) ^ 16(2p + P)4 ^ 



4P|(p) + 16P4(P) 4p2(p); il^ 4p2(p) +2P2(P) 

- 2m ) + ^^'^ ^^^^ + + ^''^ 1^ I (2pTpF " WTW ' ^ 

+ Pi«. + Pu.) - 64(^^f^Pi(Pi.Pu) - - ^(Pi. + Pi. + Pi.)Pi 



8P2(p)3 F2(p)2 y V 4P2(p)2 2F2(P) 

15 lOP? 2P,4 30 2 , -P? 3 , / 10P2 8P2 



'2F^{p) 32p3 (4p)4 16p2 P3(p)46p2 ' Ap' \F^{p) 4pF^{p) 

2 /I d2 d2 q \ / Z32 Q \ O n4 / gp2 



P2 4P2 p2 3 W p2 3 \ / 



Ap' Fi{p)H6p^ 4pVV4P-4^ '^B-2^^) Fi{p)U{B-^^f 



Is 



Pfs 15 



-I i± I- 
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32p3 I6(2p)4 16pVV 6^2 6^2(232 +2p) m 2md2 + 2p) 



6V2{^ + 2p) '6^/2 2md4^^ 7 \A{^ + 2pr 24V2{^+2p)^ 

i-Pld^-,3d^^' 90. , 2940c , 72 105p 



6^2(2^^+2^)2 48^2(2^^+2^)3 4,71(2^^+2^) m{^+2py 



TO(2^+2p)3; Vl6p' 4p;V(23^ + 2p) 6A/2(2i2+2p) 

24c 120c 70c 28p SOp^ 



6v^(2a^ + 2p) 24^2(2^+2^)2 ' 6^{^+2p)^ m{^+2p) ' + 2py 

?p2 8c\/ Pf 21P4 105 lOSPfXI TT I 

^ ~ eyiy I ~ (4pF ^ (4pF ^ (4pF " JW) r2p^'''''P^^^' 

TT 3 TT \ -1 / — P? 

2a, + 2p + P 2p + B-^^^^^^ 4(2^ + 5- (^3^4^) 

r.o , 2i?4 4P2 _5P2 ^ p4 



(2^. + 2p + P)4 {^ + 2p + By"A{2p + B- ,^fl^r,, f 16(2p + B 



15 , , 10P2 2 -P2 

+ , , + 1 , o.. , „ )(- 



^(^^+^- (^+2;+B) )^ ^(2^^+2^+5)3 2^+2p+P^^4(2p+P- ^_^f2;+B) )^ 

2i2p + B-^-^;)^^ 2(^11 + By^' -^^^ 
F2(p)= 23^+5-4(1^, and ^^3(p) = ^ + 2p + i?. 

Now from this eq.(Cll) the elements of transition matrix (T) can be found. As shown in ref. [21], the 
gluonic excitations are orthogonal to the ground states i.e. 

(1 I 1*) = (1* I 1) = (2 I 2*) = (2-^ I 2) = 0. 

This gives 

Til* = 2^1*1 = T'22* = 12*2 = 0. 

Thus the only possible transition amplitudes are Tn, T12, T21, P22, P12* and T21*. As the eq.(C5) and 
(C6) are similar (interchanging Ri and R2 in eq.(C5) gives eq.(C6)), so 

Til ~ T22, 
T12 = T21, 

and 

7l2* = 721* • 
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